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ABSTRACT 

At large magnetic Reynolds numbers, magnetic helicity evolution plays an important role in astro- 
physical large-scale dynamos. The recognition of this fact led to the development of the dynamical 
a quenching formalism, which predicts catastrophically low mean fields in open systems. Here we 
show that in oscillatory afl dynamos this formalism predicts an unphysical magnetic helicity transfer 
between scales. An alternative technique is proposed where this artifact is removed by using the 
evolution equation for the magnetic helicity of the total field in the shearing-advective gauge. In 
the traditional dynamical a quenching formalism, this can be described by an additional magnetic 
helicity flux of small-scale fields that does not appear in homogeneous a^ dynamos. In afl dynamos, 
the alternative formalism is shown to lead to larger saturation fields than what has been obtained 
in some earlier models with the traditional formalism. We have compared the predictions of the two 
formalisms to results of direct numerical simulations, finding that the alternative formulation pro- 
vides a better fit. This suggests that worries about catastrophic dynamo behavior in the limit of large 
magnetic Reynolds number are unfounded. 
Subject headings: MHD — turbulence — Sun: magnetic fields 



1. INTRODUCTION 

While the possibility, and indeed need, for astrophys - 
ical dynamos was recognized quite early ()LarmoHll919D . 
the study of dynamos has since been troubled by a 
number of pro blems. Cowling's anti-dynamo theorem 
(jCowlingl 11933) initially appeared to demonstrate that 
the entire concept was impossible, though Parker (1955) 
eventually discovered the physics behind what has come 
to be called the a effect. Cowling's anti-dynamo theorem 
was finally shown to be largely inapplicable by analyti- 
cally solvable dynamos such as the Herzenberg dynamo 
(JHcrzenbcrg 1958). Once the possibility of dynamo ac- 
tion was demonstrated, t he development of mean -field a 
dynamo theory followed (jSteenbeck et al.|[l966D . which 
describes the generation of poloidal field from toroidal 
fields. 

While the generation of toroidal magnetic fields from 
sheared poloidal fields is straightforward through the ft 
effect, the reverse process is tricky. Without it however, 
dynamo action is impossible. The a effect, which re- 
lies on helicity (twist) in the fluid motion, allows for the 
generation of strong large-scale magnetic fields such as 
those observed in the Universe. It can drive dynamo ac- 
tion on its own (a^ systems), but as shear is ubiquitous 
in astrophysics, shear- amplified dynamo action is gener- 
ally expected to outperform a^ dynamos. Accordingly, 
afl dynamos, which combine the effects, are expected to 
be the dominant type of natur al astrophysical dynamo 
(jHubbard fc BrandenburjIMl . 

More recently however, ther e were indi cations, first 
suggested by IVainshtein & Catt aned (|1992| ) , that the a 
effect decreases catastrophically already for weak mean 
fields in the limit of large magnetic Reynolds num- 
ber (i.e. low non-dimensionalized resistivities). Such 
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behavior would imply that mean-field dynamos driven 
by the a effect could not generate the observed large- 
scale magnetic fields. This claim stymied the field 
of large-scale dynamos for the 1990s. While strong 
fields are observed in nature, the theoretical under- 
standing appeared to have been cut down. Eventu- 
ally it was recognized that this behavior is not gener- 
ally applicable, being restricted to two-dimensional sys- 
tems, or to homogeneous (non-dynamo generated) mean 
fields (JBlackman fc Brandenburg] 1200^ . a nd large-scale 
dynamo simulations became com mon (jBrandenbursJ 
120011: [Brandenburg fc DobleH 120011) . These new sim- 
ulations occurred alongside the realization that mag- 
netic helicity conservation, through the dynamical a 
quenching formalism, provides an excellent theoreti- 
cal understanding of the saturation of a-effect dy- 
namos: the build-up of s mall-scale magnetic hel i- 
city quenches the a effect (jField fc Blackmaiil I2002D . 
Even so, the question of catastrophic quenching has 
remained open, with indications of saturated large- 
scale field strength decreasing with increasing mag- 
netic Reyno lds number for shearing sheets and open 
a^ systems (jBrandenburg fc Subramanianll2005"bl ). Fur- 
ther, while the saturation field strength in a^ sys- 
tems with periodic or perfectly conducting boundaries 
has been found to be independent of the resistivity 
for adequately (and in practice modestly) super-critical 
ReM, the timesc ale to reach satura tion increases lin- 
early with ReM (|Brandenburd l200ll ) . This has led to 
the study of magnetic helicit y flux e s (|Vishniac fc Choi 
12001; Brandenburg fc Sandi^ 120041 : iMitra et all 120101: 
ICande laresi ct al. 20 ll), where the hope is that, because 
the build-up of small-scale magnetic helicity quenches the 
a effect, stronger and faster growing dynamos should be 
possible if the helicity is, instead, exported (as it can- 
not be destroyed except through the action of true, i.e. 



microphysical, dissipation). 

Probing the reality of catastrophic quenching is nat- 
urally difhcult. Analytical theory is impossible, and di- 
rect numerical simulations are limited to Rejv/ that, while 
significantly super-critical for many systems, are never- 
theless orders of magnitude below those of astrophysical 
systems. The dynamical a quenching formalism allows 
probing large Rcm in systems it can handle, but its va- 
lidity there cannot, of course, be directly verified. While 
the evidence for and against catastrophic quenching is 
limited, resolving the issue is a crucial step in advancing 
dynamo theory. 

The continued improvement in techniques to measure 
turbulent dynamo coefficients from simulations has en- 
abled new approaches to evaluating different formula- 
tions of the dynamical que nching formalism. In partic- 
ular, the test-field method fSchrinn er et al.l 12005'. '2007") 
has been used to rule out the possibility of catastrophic 
quenching of the turbulent magnetic diffusivity rji in a^ 
dynamos (jBrandenburg et aLll2008aD . Recent advances 
in the th eory of magnetic helicity fluxes in the presence 
of shear (jHubbard &: Brandenburgll2011[ ) have led us to 
continue these developments in dynamical quenching by 
revisiting earlier results from shearing systems. Some- 
what surprisingly, these developments return the 1-D a 
dependent models to the first 0- D a dependent models 
(jBlackman fc Brandenburg l2002f) . In addition, we shall 
extend here earlier numerical studies of a^ dynamos in 
open systems. 



2. MEAN-FIELD MODELING 
2.1. Mean-field dynamo action 

We reproduce here some basic results of mean-field 
modeling. The dynamos we will consider are in the fam- 
ily of a^, aO, and a^fl dynamos, i.e. dynamos where 
the conversion of toroidal field to poloidal field occurs 
through the a effect, while the conversion of poloidal to 
toroidal field occurs through the a effect, the 51 effect, 
and a combination of the two. In practice, because some 
conversion of poloidal field to toroidal field through the a 
effect is always present, aQ dynamos are an approxima- 
tion in the limit that the il effect is much stronger than 
the a effect. All three dynamos, in an infinite, (shearing-) 
periodic system arise from the same eigenvalue problem. 

Although we will focus in this work on the discussion of 
results from numerical simulations, these results are bet- 
ter understood in terms of linear theory. We assume a 
standard, isotropic homogeneous a, turbulent resistivity 
rjt, and consider a system with shear velocity Us = Sxy. 
We make a standard mean-field decomposition using xy- 
planar averaging throughout, with over-barred upper- 
case variables denoting averaged quantities and lower- 
cased non-overbarred referring to fluctuating quantities, 
e.g., 



netic helicity h = A ■ B, where we use 
1 



^x^y Jx. 


/ A- Bdxdy, 
Jy 


(3) 


hjn= A- B, 




(4) 


hf = h- h„i - 


= a ■ b, 


(5) 



i.e., h is the xy- averaged magnetic helicity density and 
h„i is the magnetic helicity density carried by the large- 
scale fields. Note that while hf is the magnetic helicity 
density carried by the small-scale fields, it is still a mean 
quantity. We will further define kf as the wavenumber 
of the energy-carrying scale of the turbulence, fci as the 
scale of the mean-fields, and B^q as the equipartition 
magnetic energy while working in units for which B^q = 
Urms, the rms turbulent velocity. 

With these definitions and averaging choices, as long 
as f/ = 0, the mean field equations are written as 



(6) 

(7) 
(8) 



where J = V x i? is the mean current density in units 
for which the vacuum permeability is unity, and £ is 
the mean electromotive force, which is here expressed 
in terms of isotropic a-effect and rjt is the turbulent re- 
sistivity (jBrandenburd [200TI ). Accordingly, the mean- 
field problem for a one-dimensional a^ft dynamo with 



£ = 


aB 


-ritJ-: 


dA 
dt 


= £- 


-vJ, 
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f^x 



and kz = fci (i.e., averaging over the xy- 



plane) reduces to the eigenvalue problem 



XB 



-mkl 


— iafci 





iaki -\- S 


-mkl 











-mkl 
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(9) 



where rjT = rj + rjt is the total, microphysical and turbu- 
lent, resistivity. The growing mode has eigenvalue and 
eigenvector 



X=\aki\^l-iQ-TjTkl, (10) 

B = Bo f sin fciz, sgn (afci)(l + Q^fl^ sin(fciz -f </>), oV 

(11) 
where 

(12) 



^-4 

ak\ 



is a measure of the relative shear and 
^ — l/2arctanQ 



(13) 



B = B + b, 
B ' 



^x^y Jx Jy 



B dx dy. 



(1) 
(2) 



An important deviation from this notation is the mag- 



is the phase between B^ and By. The growth rate of the 
a^n mode is 



ReA = |aA:i|^(l + (1 + Q2)i/2)/2 - ^^kf. (14) 
From the above, we draw some significant conclusions 



true for both ail and the more general a^il fields: 



\A-B\ 



(1 + Q')'/' 2 . 

Bn smc 



fci 



(15) 



B' 



Bl fsin^ kiz + (1 + Q^)i/2 sm^[kiz + (j)]\ , (16) 

i.e., the magnetic helicity density and the current heli- 
city density Cm, = J ■ B oi the mean-field are spatially 
uniform, while the amplitude of the mean-field is not 
spatially uniform ii S y^ (Q 7^ and (j> 7^ 7r/2). 

We next make the aQ approximation, assuming that 
IQI ;» 1. We also consider only the case of a, fci , 5 > to 
simplify notation (the other cases are analogous). This 
implies that 



(1+Q2)1/4^q1/2»i^ 
^ 4' 

ReA= \akl\^/QJ2 - rjrkl 



(17) 
(18) 
(19) 



At constant 77T and S then, the system will be stationary 
for a = ac such that 



2%fel /-)l/2 _ /I ^ 



Further, in such a state we have 



K 



2r]Tki 



(B2) 



(20) 



(21) 



The a effect from maximally helical turbulence has a ~ 
Tytfc/, so the mean magnetic field of an aQ dynamo is 
expected to have very low helicity. As we will see, this is 
an important consideration. 

2.2. Catastrophic a-quenching 

Given the level of interest, it should be noted that 
"catastrophic" a-quenching has not been consistently 
defined. We will choose the following definitions: 

• Type 1 catastrophic quenching is probably the 
most extreme case. Here, the saturated mean-field 
strength varies inversely with Rcm (or some non- 
negligible negative power or similar). 

• Type 2 catastrophic quenching is well understood 
in an a^ dynamo in a triply-periodic setup as dis- 
cussed in Section 12.31 Here, the time required for 
final saturation scales linearly with Rcm (or some 
non-negligible positive power thereof). 

A well known example of Ty pe 2 catastrophic que nching 
is seen in the simulations of iBrandenburgI ()2001h . while 
Type 1 catastrophic que nching has been suspected to oc - 
cur in the simulations of | Bra ndcnbu rg fc Dobleii (|2001[ ): 
[Brandenburg fc Subramama n (20051^, but this will be 
challenged by the present work. 

It should be noted that both Type 1 and 2 quench- 
ings might be less than fully catastrophic in practice. A 
system which rapidly reaches an RcM-independent field 
strength and then resistively decays could be Type 1 
and yet have a significant field for all relevant times. 
Similarly, a system could take a prohibitive resistive 



time to fully saturate, but already reach significant field 
strengths on dynamical times. 

Given the name a-quenching, it would be appropri- 
ate to define a quenching type based on the value of 
a. Such a definition is quite difficult however, as in the 
saturated regime the dynamo-driving effect counterbal- 
ances resistive decay, so the net dynamo-driving terms, 
including the turbulent resistivity that must accompany 
an a-effect, are expected to vary with rj (and so with 
Re^/). 

2.3. Dynamical a-quenching 

Dynamical a-q uenching is a theoret i cal ad vance, first 
introduced by iKJeeorin fc Ruzmaikirj (Il982f) a nd more 
recently seen in iBlackman fc Brandenburg (2002), that 
uses the magnetic a-effect of .Pouquet et al.. (1976, ). Un- 
der that hypothesis, the actual a effect in a system can 
be decomposed into a component due to the kinetic ef- 
fect, aK, and a component due to the backreaction of 
the magnetic fields on the flow, a^/: 



a ^ aK + aM, olk 



OLM 



3p 



jb. (22) 



The mean current helicity density of the small-scale held, 
j ■ 6, is not a tractable quantity, but in general it is well 
approximated by the mean magnetic helicity density of 
the small-scale field, hf = a ■ b, through j ■ b ~ k^^a ■ b; 
see IMitra et al.l ()2010l ) for details and results in an in- 
homogeneous system. Recall that under this definition, 
small-scale magnetic helicity is a mean quantity. 

The mean small-scale magnetic helicity can be found 
by subtracting the evolution equation of the large-scale 
magnetic helicity from that of the total helicity. This 
can be determine d from the uncurled induction e qua tion 
[see Section 3 of iBrandenburg fc SubramanianI (pOO Sa). 
noting the sign error for the V0 terms in their Equations 
(3.33) and (3.44)], 



E = -U X B + 7]J, 

dA ^ „, 

B = V X A, 



(23) 
(24) 
(25) 



where E is the electric field. After some vector identities, 
we arrive at 



dh 

dhm 
dt 

dh 

dt 



-- -2r]J B V T, 
= 2£ • B - 2r; J • S - V • !^„ 



= -2£B-2rjjb-V J^f, 



(26) 
(27) 
(28) 



where JF — !Fm + ^f is the sum of large-scale and small- 
scale magnetic helicity fluxes and £ = u x b: see Eq. ([5]). 
Note that the contribution of {U x B) • S = in Equa- 
tion (1^51) is split into flnite terms of opposite sign ±£ ■ B 
in Equations (P7|) and (P5|) . The gauge term in Equa- 
tion (p4|) is in cluded in the flux terms: for a complete 
discussion see iHubbard fc Brandenburg! ()2011[ ). Using 
J • 6 ~ fc?a • b. Equation (P5| can be evolved in a mean- 
field simulation if a form for the flux term is assumed. 



We call this traditional dynamical a-quenching. In ho- 
mogeneous, periodic systems, such as homogeneous a^ 
dynamos in triply periodic cubes, the flux term vanishes, 
and the concept behind dynamical a-quenching can be 
tested. The application of dynamical a-quenching to 
this system predicts Type 2 quenching: there is an expo- 

nential growth phase which ends when B /Bl = ki/kf 
(jBlacknian fc B randenburg 20021 . Subsequently, there 
is a resistively controlled saturation phase with time 
\/2rjk\, finally ending at a saturated field strength of 

B^/Bl^ = kf/ki (IBrandenburg|l2Q0l . 

Recent work suggests that the appropriate ansatz 
for the flux of mean small-scale magnetic helicity 



is diffusive, 
( 



Hubbard fc Brandenburg 



Hubbard fc Brandenburg 



with sub-turbulent diffusion coefficients 
However, recent work 
has also demonstrated 



2010) 



2011 



that shear poses a unique problem which can be seen in 
the case of a shearing-periodic setup at a moment when 
all quantities are periodic except for the imposed shear 
flow Us — Sx y. In that case, the helicity flux has a hor- 
izontal component, (Us x B) x A, which is not periodic 
and has a flnite divergence. While the existence of this 
net flux through the shearing-periodic boundaries might 
be unexpected, the need for it can be simply explained. 
The solution of an a^Q, dynamo has spatially uniform 
large-scale helicity, as quantifled by Equation (|15l) . but 
the S ■ B term in Equation (P?]) depends on z. A flux 
term with a flnite divergence is required to balance the 
equation. This flux term follows naturally from the re- 
quirement that whatever terms the mean electromotive 
force produces in the evolution equation for the mag- 
netic helicity of the mean field, it should not affect the 
evolution of magnetic helicity of the total field. In other 
words, no terms involving £ should appear in the evolu- 
tion equation for A ■ B + a ■ b. Any term with £ in the 
equation for hm should thus be absorbed by such a term 
with opposite sign in the equation for hf. 

To elucidate this further, let us consider the equation 
for dtA — £ for the mean field. Dotting this with B gives 
the contribution 5 • B for the production oi A- B. We 
still need the contribution from A ■ dfB, i.e., A ■ V x 5. 
Using the identity 

A-Vx£ = £-VxA + V-(£xA), (29) 



we have 



^A-B = 2£-B + V-(£xA) 
at 



(30) 



where dots indicate the presence of other terms not in- 
volving £ for the full equation. Thus, the evolution equa- 
tion for a ■ b must then be of the form 

d- 



dt 



ab = -2£ -B-V ■{£ X A) 



(31) 



so that the evolution oi A ■ B + a ■ b is not effected by 
the £ terms. In the traditional dynamical a-quenching 
formalism, this was only true of the 2£ ■ B term, but the 
divergence of 5 x A had been ignored. 
Allowing now for all the other terms in Equations (|27p 



and ((28| , our full set of equations is 
dh„ 



dt 



2£-B-27jJ-B-V ■{T„,-£xA), (32) 



^^-2£.B-2ru-b-V-{J^-J^'^ + £xA), 

(33) 

where .^„ is the resistive component of ^m . When B 
takes the form in pT|) and £ = aB — r]tJ, the V • (5 x A) 
terms cancel the £B terms. IBerger fc RuzmaikinI ()2000f ) 
have estimated that this £ x A flux can be important in 
the Sun. 

If the flux term is not correctly handled, we can expect 
the generation of artificial helicity "hot-spots" through 
the £ ■ B terms, which will nonlinear ly back- react on the 
dynamo through Equation ([22l) . While an adequate dif- 
fusive flux may be able to smooth out such, this poses 
a clear potential difficulty in applying dynamical a- 
quenching to shearing systems. 

If a mean-field model is solved in terms of the mean 
magnetic vector potential A however, then hm is known 
at every time step. Thus, rather than evolving Equa- 
tion ([28l) . one can evolve Equation (f26|) to find hf — 

h — hm, avoiding the £ ■ B terms. One known difficulty 
with this alternate technique is that spatially homoge- 
neous components of A may develop and cause spurious 
spatial variation in aM when the latter is defined in terms 
of A • B. This homogeneous component arises from nu- 
merical noise: as a constant A is curl-free, physically 
motivated equations cannot generate or, unfortunately, 
erase it. Accordingly, we must artificially treat the issue 
by subtracting out the volume averaged {A)v- We refer 
to this technique of calculating hf as alternate dynamical 
a-quenching. For systems with no native spatial varia- 
tions in a, nor any instabilities in the spatial variation of 
a, this procedure will in practice return one to the first 
attempts to apply dynamical a-quenching using volume 
averages (jBlackman fc Brandenbur^l2002[ ). 

We use an a^ dynamo to test alternate dynamical 
a-quenching against traditional dynamical a-quenching 
(which, in this system, should be identical as there are no 
spatial variations and so no fiuxes). We show the agree- 
ment in Figure [T] The small difference that develops is 
due to a smaller rms spatial noise of aM in the alternate 
quenching case. 

2.4. Investigation procedure 

Catastrophic a-quenching lives in the high Rcm 
regime, beyond the reach of current direct numerical 
simulation or laboratory experiment. This makes con- 
firming or disproving its existence impossible. The ev- 
idence for its exis tence lies largely on mean-field sim- 
ulations (see, e . g.. iBrandenburg fc SubramanianI [2005bl : 
iGuerrero et al.l l2010f ). which confirm Tvpe 2 quench- 
ing for homogeneous isotropic periodic a^ dynamos. 
Further, mean-field simulations using traditional a- 
quenching have strongly suggested the existence of Type 
1 quenching for shearing systems. 

While we cannot simulate a^il dynamos at high Kcm, 
we are in a position to run modest Rcm shearing simula- 
tions to compare the predictions of traditional quenching 
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Fig. 1. — Mean-field simulations for an a^ dynamo at Rejv/ = 
10'^, comparing traditional (yellow/gray, solid/thick) and alternate 
(black/dashed) dynamical quenching models in a system where 
they are formally identical. 



(with and without diffusive magnetic hehcity fluxes) with 
those of alternate quenching. In the latter case, we do 
not include uncertain diffusive fluxes because the mag- 
netic helicity and therefore a-effect are not expected to 
exhibit spatial dependencies, which we confirm. 

Our procedure then is to run a direct numerical simu- 
lation of an a^Q. dynamo, extract the spatial dependency 
of a and compare it with the results of mean-field theo- 
ries. Once mean-field theories have been weighed against 
the evidence, we move to large ReM and examine the ev- 
idence for or against Type 1 and 2 quenching. 

3. NUMERICS 

We perform mean-field numerical simulations for a 
shearing sheet, with Us = Sxy, and averaging per- 
formed over the xy plane, so mean quantities are only a 
function of z, reducing the problem to a one-dimensional 
one. Our mean field equations are evolved using the same 
algorithm as the Pencil Code (see below), but due to 
ease of implementation at the time, and low numerical 
load, run using the Interactive Data Language (IDL). 

We formulate the mean £ through the standard for- 
mula (Et, where rtt is assu med not to be quenched; see 
iBrandenburg et aLl (j2008a| ) for a numerical justification 
of this. The total a is given by the sum of the kinetic 
aK, presumed constant, and the magnetic aM- Accord- 
ingly, da/dt — daM/dt. We solve the two systems of 
equations 



da 

dB 

'dt 



-2mk^ 



£ ■ B a — ax 

1 — 

B^q Vt/V 



V X {£- T]J) 



VaV^a, (34) 



(35) 



for traditional quenching, see, e. g.. Equa t ions (9 .14) and 
(9.15) of Bran denburg fc SubramanianI (|2005aD . where 
the helicity fluxes have been c a,st in the form of diffusion 
terms following the results of iHubbard fc Brandenbur"^ 
120101 where it was found that the flux was proportional 
to the gradient of the magnetic helicities. The diffusive 
helicity flux has diffusion coefficient 2?q which will be 
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Fig. 2. — Top left panel: By in a frame comoving with the 
dynamo wave. Top right panel: a^x- Middle panels: butterfly 
diagrams of components of o(z,t). Bottom panels: Time-averages 
of the middle panels. Note differing y axis scales, whic h im plies 
that the quenching of ayy is nearly uniform. See Section |4. II 

scaled to rjt- Alternate quenching solves instead 



dh 

at ~ 


-2rj ( J • B + aMB^^jTjt) , 


(36) 


aM = 


--i^tkj{h-A-B)lBl^, 


(37) 


dA 

at 


^I-TJJ. 


(38) 



Note that for alternate quenching we also enforce 
J Adz = at every timestep to avoid drifts in the 
magnetic vector potential. The essential difference be- 
tween the two approaches can be traced back to mutually 
canceling contributions to the large-scale and small-scale 
magnetic helicity flux of the form ^£ x A. 

Our direct numerical simulations are made using the 
Pencil Code, a finite-difference scheme sixth order in 
space and third order in time. In the Pencil Code 
runs, we use the test-field method (TFM) to determine 
components of the a ten sor as a function of position . For 
information on TFM , sec Brande nburg et al.l (|2008b[ ) and 
iRheinhardt fc Brandenburg (2010ll 

4. MEASURED a PROFILES 

4.1. Direct simulation 

In Figure [5] we present data for the z dependence 
of a well into the saturated regime for a direct sim- 
ulation with Rcm = 27 and fc//fci — 3. As in 
IHubbard &: Brandenburg! (|2011| ) , the butterfly diagrams 
are shifted to the frame comoving with the traveling dy- 
namo wave, as demonstrated in the top-left panel. This 
allows us to take meaningful time-averages while retain- 
ing spatial information. In the top-right panel we show 
the volume rms of a^x in a semi-logarithmic plot, which 
demonstrates that the system (including the small-scale 
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Fig. 3. — Butterfly plots of magnetic helicity. Top panel: 
{A)v ■ B, i.e., the fictitious component of hm due to a spatially 
homogeneous component of .A. Middle panel: /i^, adjusted for the 
top panel. Bottom panel: hf. While there may be some spatial 
structure in the bottom panels, it is intermittent in time, and the 
residual from a near-cancellation (the bottom two panels use a very 
different scale than the top one, see the color bars). 

fields in the TFM) is in a steady state for the time 
interval consid ered. The deep spike marks a reset of 
the t est-fields (|Ossendriiver et aLl 120021 : iHubbard et al.l 
I2009D . The middle two contour plots show the two im- 
portant components of a in the comoving frame. The 
middle left panel shows axx which aids the J7 effect in 
converting the poloidal Bx into By, and the middle right 
panel the vital ayy which provides the conversion of the 

toroidal By into Bx- 

The bottom panels are time-averages of the middle 
panels. It appears from the contour plot that ayy shows 
spatial variation, which is confirmed when a time aver- 
age (in the shifted domain) is taken as seen in the bottom 
right panel. However, in the bottom left panel it is clear 
that the actual result is that ayy is strongly quenched 
compared to axx- This implies that the spatial variation 
seen in ayy is merely spatial variation in the residual a 
effect: the quenching itself is nearly uniform. 

Figure 131 for the same simulation, shows the difficulty 
mentioned in Section 12.31 namely that a spatially ho- 
mogenous component of {A)v can generate a spurious 
magnetic helicity signal. We must also note here that 
the quenching is blatantly non-isotropic. A study of this 
effect is beyond the scope of this paper: we expect it to 
be a full project in its own right, and intend to study it 
as such. 
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Fig. 4. — Field and a = an + dim (taken at the final time) 
profiles for mean-field simulations, approximately saturated. The 
simulations are intended to be compared with that of Figure [2] 
so it has the same Rejv/ = 27 and kj/ki = 3. Black/solid: alter- 
nate quenching formula. Blue/dashed, red, black and yellow/dash- 
dotted use the conventional quenching formula with "Dajm = 
10, 1, 0.1, 0.01, 0, respectively. 

4.2. Mean-field approaches 

In Figurc|4]we show energies and aM profiles for mean- 
field simulations similar to that of Figure [5] (Rcm = 27, 
kf/ki — 3). The mean-field simulations use traditional 
quenching with T^a/flt ranging from to 10 and a run 
with alternate quenching. None of the traditional models 
match the uniform quenching that is measured in Fig- 
ure [21 showing large spatial variability that derives from 
the £ ■ B term in Equation ([28]) . not even the model with 
Da = lO^t- The decrease in spatial variation of a with 
increasing "Da suggests that the traditional model could 
be made to function with an adequate diffusion term, but 
this term would need to be absurd in scale (and would 
hopelessly distort any simulation with "real" spatial vari- 
ation in aM that needs to be correctly captured). The 
alternate quenching formalism does result in the uniform 
quenching, which is unsurprising as it eliminates the spa- 
tial forcing from £ ■ B- 

We take this as strong evidence that the alternate 
quenching formalism is superior to traditional quench- 
ing in sheared systems where drifts in A are tractable - 
and that results obtained with traditional quenching in 
the presence of shear should be viewed with suspicion. 

5. MEAN FIELD: LARGE MAGNETIC REYNOLDS 
NUMBERS 

5.1. Early times 

For early times, the predictions of both dynamical a- 
quenching formalisms predict behavior similar to that 
of a^ dynamos: exponential growth of the mean fields 
(and corresponding growth of au) until the total a effect 
is reduced enough that the growth rate is reduced to 
a fraction of its original self. This occurs when |a| = 



\2ri^kl/S\, i.e., when 



^o--^) 



\aK\ - |27?4fei/^l 



(39) 
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Fig. 5. — Fie ld magnitudes for mean-field simulations discussed 
in Section 15.11 into the non-kinematic regime. The simulations 
use Re A/ = 10'^, kf/k\ = 3. Black/solid: alternate quenching 
formula. Blue/dashed, red, black and yellow/dash-dotted use the 
conventional quenching formula with T>a/Tit = 10,1,0.1,0.01,0, 
respectively. Note the oscillations, which persist at some level even 
with Va = 0.1»;t. 

In terms of magnetic hclicity, this becomes 



{a-b)\ 



\aK\ 



\2r,?,kl/S\) . (40) 



Using th e standard appro ximations for fuhy hehcal tur- 
bulence ()Sur et al.ll2008( ). namely t ~ 1/urmskf, ax — 



Urms/3 and % 
reduces to 



g/3, and writing B, 



cq 



s, this 



{a-b)\ 



2klB^. 



-Scq 



iB\S\ / kf 



-^. (41) 



As the growth is rapid, we will have ft,,, 
this stage, and so 



1 



2,kj\S\ 



Bog 



hf during 



(42) 



However, afl dynamo mean-fields are only weakly helical, 
i.e. hm <C B'^/ki. Under the assumptions that the mean 
field is approximately stationary, and that the shear is 
strong enough to use Equation (|2ip as an approximation. 
Equation (|42p implies that: 



(W 



2r/Tfc' 



'^7n '^ 



2aKki 




(43) 



As we have made the afi approximation that \S\ ^ 
jaiffcil, this implies that an aO, field first feels nonlinear 
effects for mean-field energies that are already in super- 
equipartition. 

In Figure [5] we show the early evolution of a mean- 
field dynamo with Rcm = lO'^, uk = —1/3, S — 1 and 
kf/ki = 3. Equation ()43p implies that the exit from ex- 
ponential growth occurs for B w l.lSScq, which is well 
captured by alternate quenching, and traditional quench- 
ing with strong diffusive fiuxes. 



5.2. Late times 
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Fig. 6. — Field magnitudes, for the same mean-field simu- 
lations as Figure [5] {Rejj/ = 10 , kf/ki = 3), into the non- 
kinematic regime. Red/solid and over-plotted black/dashed: al- 
ternate quenching formula for Rejv/ = 1000 and 2000 respec- 
tively. Blue/dashed, red/dash-double-dotted, black/dash-dotted 
and yellow/dash-double-dotted use the conventional quenching for- 
mula with 'Da/'Tjt = 10,1,0.1,0, respectively. The curve for 
'Da/vt = 0.1 is extremely strongly coarse-grained for visibility: 
the oscillations (for that run) are in fact far more frequent than 
shown and would be a solid band if plotted fully. 

We can analytically estimate the final field strength 
of the dynamo for the alternate quenching formalism, 
while for traditional models the problem is nonlinear as 
can be seen in Figure |4l The final state is achieved when 
dh/dt = 0, i.e., when J ■ B = {j ■ b). Combining this 
with Equation (j2ip and assuming that the shear is strong 
enough that a must be fully quenched, a = uk^-olm — 0, 
we find 



b' 



\S/2aKki\{kf/k^fBl 



(44) 



In Figure [6] we show the late time evolution of the same 
mean field dynamos as in Section 15.11 It is clear that, 
without significant {T>a > O.lr^j) hclicity diffusion, the 
solution for traditional quenching is unstable and drops 
to resistively small values. This is not surprising as 
the problem becomes highly nonlinear. However, with 
moderate diffusion the field strength behaves smoothly, 
with the final energy level increasing with diffusion co- 
efficient. Even so, the saturation level of the traditional 
quenching model with Da = lO??* is significantly below 
that of the alternate quenching model. While the dif- 
fusion does smooth out the helicity hot-spots, the spa- 
tial fiuctuations of aM in Figure [3] have a noteworthy 
impact on the final dynamo state. Finally, the satura- 
tion level of the alternate model matches the estimate 
from Equation (pi|) of B ~ 3.7i?cq, to within the limit 
that an adequate residual a is needed to sustain the 
field against turbulent resistive decay. Further, the over- 
plotted black/dashed alternate-quenching curve is for 
Rcm = 2000, double that of the red/solid alternate- 
quenching curve. The overlay implies that we have 
reached an asymptotic state independent of ReM , which 
is in agreement w ith earlier work assuming perfect spa- 
tial homogeneity (jBlackman fc B randenburg' '2002^ and 
with simulations (see Figure 6 of lKapvla fc Brandenburg 
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Fig. 7. — Time series for a^ dynamos in open systems. 

Black/solid: potential field extrapolation in the vertical direction. 

Red/blue/dashed: vertical field condition on the vertical direction. 

Potential field and VF5 have Re^^f = 156 while VFc has Hbm = 86. 

1 /2 
The dash-double-dotted line corresponds to 1/kJ , i.e. the energy 

level associated with the end of the kinematic phase. 
[2OO9I) . 

6. DIRECT SIMULATIONS OF OPEN SYSTEMS 

Numerical resources limit our ability to probe the high 
ReM regime. However, we have run three simulations of 
a^ dynamos in an open system, i.e., a system which can 

export magnetic helicit y. ^This system is the same one 

as considered in Brandenburg & Subramanian (2005b); 
a helically forced cube, periodic in the horizontal direc- 
tions and with vertical field conditions in the vertical 
directions, which we have run for Rcm = 86 and 156. 
Additionally, as the vertical field condition is frequently 
used instead of a proper vacuum condition, we also per- 
formed a Rca/ = 156 run with potential field condition 
in the vertical directions. The resulting time series are 
given in Figure [T] Our resolution was 128'^, for runs 
with Winax — 0.15, M,ms — 0.05 and 77 = 2 x 10"'* (for 
Rejvf = 86) or ?7 = 10""' (for the other two). The veloc- 
ity boundary has a stress-free vertical condition, and the 
entropy a symmetric one. 

Unlike the results ^reported in 

iBrandenburg fc: SubramanianI ()2005b[ ). there is no 
clear indication of a reduction in the strength of the 
mean field for higher magnetic Reynolds number, 
even though the runs were followed for resistive times. 
However, the use of vertical field conditions as a proxy 
for vacuum conditions appears to be a poor one. Note 
that there does not appear to be a slow resistive phase. 
This lack is expected as the open boundaries allow 
the system to export total magnetic helicity (not just 



helicity of the small-scaled field). Thus, the system 
should reach a steady state where exchanges of helicity 
through the boundary balance preferential destruction 
of small-scale helicity on dynamical times, and for small 
total helicities. 

7. DISCUSSION AND CONCLUSIONS 

We have used the test-field method to examine the 
predictions of catastrophic a-quenching resulting from 
dynamically-quenched mean-field models in shearing 
systems. Formulations for dynamical a-quenching which 
are superior for the problem of shearing systems do not 
predict Type 1 catastrophic quenching (reduced field 
strength) but do predict Type 2 quenching (long final 
saturation times) , extending results that do not allow for 
spatial variations of a (Blackman & Brandenburg 200^ 
to models that do. We have further revisited simulations 
of a^ dynamos in open systems and, at admittedly quite 
modest RcAf , found no evidence of field strength scaling 
inversely with Rcm- 

The picture we see now for a-effect dynamos, moti- 
vated by the concepts and formalism of dynamical a- 
quenching, is one of exponential growth during a rapid 
initial saturation phase. This phase ends when the mag- 
netic helicity in the small-scale fields is comparable to 
the helicity in the forcing that generates the a-effect. At 
this point, the total magnetic helicity in the system has 
not changed from its initial value. If the system is open, 
exchanges with the exterior (Section [6|) will tend to keep 
the total magnetic helicity roughly constant, and the sys- 
tem will then not evolve resistively. On the other hand, 
if the system is closed the preferential resistive destruc- 
tion of magnetic helicity of the small-scale field allows a 
further resistive growth phase. 

It is important to note that the energy in the large- 
scale field is bounded below by its helicity. Weakly heli- 
cal large-scale fields are possible, which can have super- 
equipartition fields even at the end of the kinematic 
growth phase. Weakly helical large-scale fields are a 
natural product of sheared system, so rapid growth to 
sub-equi-, equi- and super-equipartition fields are all ex- 
pected to occur in nature, although all equi- and super- 
equipartition fields in the high Rcm systems of astro- 
physics are expected to be weakly helical. 

This work was supported in part by the European 
Research Council under the AstroDyn Research Project 
227952. Alexander Hubbard acknowledges the additional 
support of a fellowship from the Alexander von Hum- 
boldt Foundation. The computations have been carried 
out at the National Supercomputer Centre in Linkoping 
and the Center for Parallel Computers at the Royal In- 
stitute of Technology in Sweden. 



REFERENCES 



Berger, M. A., & Ruzmaikin, A. 2000, J. Geophys. Res., 105, 

10481 
Blackman, E. C, 
Brandenburg, A 
Brandenburg, A 
Brandenburg, A 
Brandenburg, A 
Brandenburg, A 

400 



& Brandenburg, A. 2002, ApJ, 579, 359 
2001, ApJ, 550, 824 
& Dobler, W. 2001, A&A, 369, 329 
& Sandin, C. 2004, A&A, 427, 13 
& Subramanian, K. 2005a, Phys. Rep., 417, 1 
& Subramanian, K. 2005b, Astron. Nachr., 326, 



Brandenburg, A., Radler, K.-H., Rheinhardt, M., & 

Subramanian, K. 2008a, ApJ, 687, L49 
Brandenburg, A., Radler, K.-H., & Schrinner, M. 2008b, A&A, 

482, 739 
Candelaresi, S., Hubbard, A., Brandenburg, A., & Mitra, D. 2011, 

Physics of Plasmas, 18, 012903 
Cowling, T. G. 1933, MNRAS, 94, 39 
Field, G. B., & Blackman, E. G. 2002, ApJ, 572, 685 



Guerrero, G., Chattcrjcc, P., & Brandenburg, A. 2010, MNRAS, 

409, 1619 
Herzenberg, A. 1958, Phil. Trans. R. Soc. A, 250, 543 
Hubbard, A., & Brandenburg, A. 2010, Geophys. Astrophys. 

Fluid Dyn., 104, 577 
Hubbard, A., & Brandenburg, A. 2011, ApJ, 727, 11 
Hubbard, A., Del Sordo, F., Kapyla, P. J., & Brandenburg, A. 

2009, MNRAS, 398, 1891 
Kapyla, P. J., & Brandenburg, A. 2009, ApJ, 699, 1059 
Kleeorin, N. I., & Ruzmaikin, A. A. 1982, 

Magnetohydrodynaniics, 18, 116 
Larmor, J. 1919, Rep. Brit. Assoc. Adv. Sci., 159 
Mitra, D., Candelaresi, S., Chatterjee, P., Tavakol, R., & 

Brandenburg, A. 2010, Astron. Nachr., 331, 130 
Ossendrijver, M., Stix, M., Brandenburg, A., & Riidiger, G. 2002, 

A&A, 394, 735 
Parker, E. N. 1955, ApJ, 122, 293 



Pouquet, A., Frisch, U., & Leorat, J. 1976, Journal of Fluid 

Mechanics, 77, 321 
Rheinhardt, M., & Brandenburg, A. 2010, a, 520, A28 
Schrinner, M., Radler, K.-H., Schmitt, D., Rheinhardt, M., 

Christensen, U. 2005, Astron. Nachr., 326, 245 
Schrinner, M., Radler, K.-H., Schmitt, D., Rheinhardt, M., 

Christensen, U. R. 2007, Geophys. Astrophys. Fluid Dyn., 101, 

81 
Steenbeck, M., Krause, F., Radler, K.-H. 1966, Zeitschr. 

Naturforsch. A, 21, 369 
Sur, S., Brandenburg, A., & Subramanian, K. 2008, MNRAS, 

385, L15 
Vainshtein, S. I., & Cattaneo, F. 1992, ApJ, 393, 165 
Vishniac, E. T., & Cho, J. 2001, ApJ, 550, 752 



